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The incoherent transition strength of the exotic — e~ conversion in the ^"^Pb nucleus is 
investigated by utilizing the Continuum RPA method, appropriate for the evaluation of the rate 
that goes to the continuum of the nuclear spectrum. We find that the contribution of resonances 
lying high in the continuum is not negligible. Special attention is paid to the detailed study of 
the pronounced 1~ contribution which, according to previous calculations, dominates the overall 
incoherent rate in about all the nuclear targets. The spurious center of mass admixture to the 
partial rate originating from the 1~ excitations is explored, and its elimination is performed by 
correcting properly the dipole operators. The results found this way show that the greatest portion 
of the total 1~ contribution to the incoherent rate is spurious. 

PACS numbers: 23.40.Bw, 23.40.-s, 14.60.Pq, 21.60.Cs 



I. INTRODUCTION 

The exotic neutrinoless conversion of a bound muon to 
an electron, 



+ [A, z) ^ e- + [A, zy 



(1) 



is an interesting lepton flavor violating process 0,0,0,0. 
One of its basic characteristics is the possibility of the 
coherent channel, i.e. the ground state to ground state 
transition 0| . Experimentally, only the ratio of the 
coherent rate divided by the total muon capture rate, 
which is the dominant branching ratio exhausting a great 
part of the total (/i^, e^) rate could be measurable; by a 
judicious choice of the target nucleus, this channel could 
be free from the reaction induced background 0, 01 ■ 

The incoherent (/i~,e~) rate is a less significant por- 
tion of the total rate and much harder to calculate, but 
its knowledge is important for determining the fraction 
of the coherent process to the total (//~, e~) rate, which 
experimentally is also an interesting quantity. The theo- 
retical calculation of the total rate requires reliable coher- 
ent and incoherent nuclear matrix elements 0, ^3 ■ 
Transitions of the reaction have been previously stud- 
ied by employing various methods, such as: (i) Closure 
approximation within the shell model and the quasi- 
particle RPA (QRPA) for calculating the average con- 
tribution of the transitions to all excited states of the 
target nucleus, (ii) a Fermi gas method utilizing a rel- 
ativistic Lindhard function to compute the sum of all 
partial rates of the incoherent channel UM, (iii) state-by- 
state calculations by using shell model 0], and various 
QRPA versions JJiJ , to construct explicitly the final 
nuclear states. 
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Within method (ii), the incoherent rate is calculated 
by integrating over a continuum of excited states of a lo- 
cal Fermi sea. Therefore, this method is not appropriate 
for individual calculations of each accessible channel, but 
it offers the advantage of taking into consideration the 
part of the rate that goes to the continuum, which is not 
explicitly included in state-by-state calculations. This is 
one of the reasons why the incoherent matrix elements 
obtained with shell mod el B and the various versions 
of RPA calculations |10 [lOj, appear to be smaller com- 
pared to those of the method (ii). On the other hand, the 
common RPA and the various refinements of the QRPA 
offer a relatively simple and detailed state-by-state calcu- 
lation of all the individual low-lying excitations induced 
by the /i~ — e~ conversion operators 0,I3- 

An important conclusion of the state-by-state calcu- 
lations is that the contribution of the 1~ states to the 
incoherent (/i~,e~) rate is very large (for most of the 
isotopes it is the maximum one) compared to that of the 
other multipolarities 0^3- The portion of the 1~ con- 
tribution was found to be about 50% for all mechanisms 
leading to the /x~ — e~ conversion. Therefore, it is essen- 
tial to properly remove possible spurious contaminations 
when describing this process and other similar ones. 

The methods employed so far for the removal of the 
spurious center-of-mass (CM) admixture from the 1~ 
contribution 0, can be classified in two categories: 
(i) Those which remove from the contaminated Hamilto- 
nian the spurious terms, i.e., those containing the CM 
position (R) and the total momentum (P) operators 
and their couplings with the intrinsic Hamiltonian Hi^t 
[T^ IT^ . The diagonalization of Hint obviously gives the 
real spectrum of the studied nucleus. In this way, the 
eigenstates of the system separate into the intrinsic nu- 
clear spectrum and the pure CM excitation which can 
be omitted. Such methods, however, are tedious and not 
practical, (ii) Those which construct first a set of pu- 
rified wave functions to be used for the diagonalization 
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of the contaminated Hamiltonian. The orthonormaliza- 
tion, however, usually necessary in these methods, hin- 
ders their use. We should also mention that a recent 
method by Bes and Civitarese which removes ex- 
actly the spurious contaminations of the dipole operator, 
shows that the spurious portion is much bigger than pre- 
viously thought and other methods give. 

In QRPA calculations, an effective elimination of spu- 
rious components from the 1~ states may be achieved by 
adjusting the parameters which scale the effective inter- 
action, so that the energy of the first 1^ state becomes 
equal to zero (equal to the purely spurious center of mass 
eigenstate). In most of the cases, this requires unphysical 
values for the parameters which renormalize the individ- 
ual particle-particle {gpp) and particle-hole (gph) chan- 
nels (usually gpp ^ 1.3 — 1.5, gph ~ 0.2 — 04). Even when 
the spurious state occurs very close to zero energy, spu- 
rious admixtures remaining at higher energies cannot be 
avoided completely. 

In Ref . ^3 Etn approximate removal of the spurious 1 ~ 
components was performed by constructing the properly 
normalized purely spurious state \S) and evaluating its 
overlap with all 1~ states involved in the chosen model 
space. This showed that mostly the lowest lying 1~ state 
is affected by the translational invariance breaking caused 
by the use of empirical single particle energies and a trun- 
cated model space in RPA. This state was considered as 
fully spurious and the others were treated as real nuclear 
excitations. After removing the spurious contributions, a 
renormalization of the interaction was required for repro- 
ducing the energy spectrum of the nucleus with the use 
of realistic two-body forces. The above method is simple 
and easy to apply, but it is not as exact as that of Refs. 

Elll. 

The purpose of the present work is to study in detail 
the incoherent rate of the fi^ — e~ conversion, by us- 
ing the Continuum RPA method. We evaluate the con- 
tribution of high-lying continuum excitations. Such an 
explicit calculation has not yet been addressed. As an 
application, we study extensively the incoherent rate for 
the heavy nuclear target ^ospj^ j^y using Skyrme interac- 
tions. In order to eliminate spurious CM contaminations 
in the dominant 1^ channel, we obtain the corresponding 
energy distributions by using properly corrected dipole 
operators, which induce the ^~ — e~ conversion intrinsic 
1~ excitations. 

The paper is organized as follows. In Sec. ^ we de- 
scribe briefly the /i~ — e~ conversion operators and the 
formalism of the Continuum RPA method. In Sec. lIIII we 
calculate the strength distributions as functions of the ex- 
citation energy. The elimination of the spurious center 
of mass contaminations is also discussed. In Sec. II VI we 
summarize our main conclusions. 



II. DEFINITIONS AND BRIEF DESCRIPTION 
OF THE METHOD OF CALCULATION 

The inclusive {fi^,e^) rate is evaluated by summing 
the partial contribution of all final states |/). For spher- 
ical or nearly spherical nuclei, the vector contribution is 
given by 

Sa^Y.( — y\(f\Oa{cifm\^ (2) 

where Oa{qf) represents the /i^ — vector- type transi- 
tion operator resulting in the context of a given mecha- 
nism mediated by a photon (a = 7), a W-boson (a — W) 
or a .Z-particle exchange (a = Z). Here qj, with magni- 
tude Qf = TTifj^ — Cfc — Ef, is the momentum transferred to 
the nucleus. Ef is the energy of the final state |/) with 
respect to the ground state |0), is the binding energy 
of the muon and its mass. The transition operators 
have the form 

A 

Oa{q) - gvfv 6ca(r,)e-"^■'■^ c^r,) = ^ + ll3aT„ 

(3) 

where Tj is the 3rd component of the jth particle's 
isospin. The parameter fv = 1.0 represents the vector 
static nucleon form factor and the normalization coeffi- 
cient gv takes the value 1/6 for the photonic case and 1/2 
for the non-photonic W boson and SUSY Z exchange . 
The value of (3a depends on the model assumed. (We have 
taken the relevant values from Ref. @.) Thus, protons 
(neutrons) contribute to a given process with a "charge" 
whose value is determined by Ca(l/2) = 1/2 + (3a/(> 
(cq(— 1/2) = 1/2 — /3a/6). In the photon and Z case, 
the isoscalar and isovector components of the transition 
operator are (almost) equally important, whereas Ow is 
predominantly isoscalar. 

By assuming that the initial and final states are of def- 
inite spin and parity, a multipole decomposition of the 
operators of Eq. Q into operators TaLM of orbital angu- 
lar momentum rank L can be carried out. For spherical 
nuclei we can assume, without loss of generality, q — z. 
Then, only terms with M = survive, for which we ob- 
tain 

TaL{q) = TaLo{q) = 5v/yV4^(2i + l) 

A 
J = i 

A phase factor (— i)^ has been omitted. The contribution 
of each multipolarity to the transition rate Sa reads 

We now rewrite the rate SaL as the integral of a suit- 
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able distribution over excitation energy: 

SaL = J dERaUE) (6) 



with 

RadE) 



y (i + 

X |(/|T„i(m^ - eb - ^/)|0) \^S{E - Ef) 



1 - 



£6 



tE' 



R'aLiE). (7) 



In the above expressions, 

R'aLiE) = ^ |(/|T,i(™^ - e, - Ef)\0)\'5iE - Ef) (8) 
/ 

stands for the "strength distribution" corresponding to 
the operator TaL{q), with q ~ — eb — E. The total 
rate of Eq. © is subsequently written as 



SaL 



£6 



Mo 



-(m^ - eb)Mi + 



1 



-Mo 



where Mk = J R'ai^{E)E^dE is the fc-moment of 
R'ahiE)- The final states |/), excited by the single- 
particle operator TaL, are of particle- hole {ph) type. 
Then, the distribution R'^j^{E), and from it RaL{E) 
and SaL, can be calculated following the standard RPA 
method. 

We consider p/i excitations, built on top of the Hartree- 
Fock (HF) ground state of a closed-shell nucleus and sub- 
jected to the ph residual interaction (HF-I-RPA method). 
In this work, the quantities introduced above are calcu- 
lated using a self-consistent Skyrme-Hartree-Fock (SHF) 
plus Continuum- RPA (CRPA) model. The HF equations 
describing the ground state are derived variationally from 
the Skyrme energy functional and solved numerically us- 
ing the code of P.-G. Reinhard 0. The CRPA is formu- 
lated in coordinate space, as described, e.g., in [T^ 
and outlined below. 

The main ingredient of the model is the ph Green 
function Gl{E) in coordinate space. In particular, we 
are interested in the radial part GujirT^r' t' ] E), which 
describes the propagation of a fluctuation (or ph state) 
of multipolarity L and energy i?, excited by the oper- 
ator Vli at the point r, t and decaying via the opera- 
tor Vlj at the point r'r' - where r or r' corresponds to 
the isospin character (proton or neutron) of the fluctua- 
tion. ViiQ) stands for one of the rank-L operators Yl, 

[yLO(V2+V'')]L, [Yl±MV-V')]l, [1l±i<8)(V+V')]l 
present in the Skyrme interaction. In practice, the radial 
coordinates are replaced by points on a discretized mesh. 
A radial step Ar and a maximum value rmax (larger than 



the nuclear radius by a factor of about 3 or more, usually) 
are introduced. Then the Green function Gl{E) can be 
represented as a super-matrix in coordinate space, isospin 
character and operator indices i. 

The operator TaL{q), Eq. Q), is a multipole operator 
of the type Yl - let us label it as Vli- Therefore, within 
our RPA model, the distribution R'^^E) of Eq. JHJl is 
given in terms of the RPA Green function G^^^{E) by 

R'aLiE) = ^^TT[Tldq)G^^,^{E)TaLiq)] (9) 

in a matrix notation. In practice, this reads 

RUiE) = 144g2,/2(2L + l)^c„(r)c,(T') 



xlm j jLiqr)Gl\tirryT'-E) 

xiLiqr')drAr' , (10) 

q — rrii^ — eb—E. The integrations are to be understood as 
numerical ones, carried out by summing over the radial 
mesh points. 

The RPA Green function is given by the equation 

Gf^{E) - [1 + GLiE)'V,,,]-'Gl{E), (11) 

which is solved as a matrix equation in coordinate space, 
isospin character and operators Vli- The ph residual in- 
teraction 14cs is zero-range, of the Skyrme type, derived 
self-consistently from the Skyrme-HF energy functional 
[T^ l20l |. In this work, spin-dependent terms and the 
Coulomb interaction are omitted from Vrcs- The radial 
part of the unperturbed ph Green function of multipo- 
larity L is formally given by: 

Gl^^irryr';E) = j^mmmMllL 



ph 



Sph - E 



^ {h\VL,\p):'r{h\VL^\p)rr ] , . 

Sph + E /■ ^ ^ 

The sign of the second term depends on the symmetry 
properties of the operators Vli and Vlj under parity and 
time-reversal transformations. With h (p) we denote the 
quantum numbers of the HF hole (particle) state and 
£ph = £p^ £h is the energy of the unperturbed ph excita- 
tion. A small but finite ImE = r/2 ensures that bound 
transitions acquire a finite width. 

The particle continuum is fully taken into account, as 
follows T^: The summation over the particle states p 
in Eq. H12|) is replaced by the summation over all single- 
particle states k. The additional hole terms in the first 
term on the r.h.s will cancel the ones in the second term. 
Next, Ek is replaced by the single-particle Hamiltonian. 
Finally, the completeness of the k states and the prop- 
erties of the particle Green function are used to replace 
the sum over /c by a closed expression. Therefore, the 
only truncation introduced is the one of the radial co- 
ordinate, r < Tniax- The latter is very well justified. 
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because the amplitude of the radial wavefunctions of the 
hole states entering expression (|12f) vanishes at distances 
much larger than the nuclear radius. 

Results derived within this model for the nucleus ^"^Pb 
are presented and discussed in the next section. 



III. RESULTS AND DISCUSSION 

In the following we will present results obtained using 
the SkM*[2H parametrization of the Skyrme force. It 
describes satisfactorily giant resonances of stable nuclei, 
and therefore it is suitable for the present study. We 
have verified that our conclusions do not change when 
the parametrization SGII [TtL is used. In order to 
test the sensitivity of our results on the interaction used, 
we have also employed MSk7 , which has a large effec- 
tive mass, thereby shifting most excited states to lower 
energies compared to the more reliable SkM* and SGII. 

For the nucleus ^"*Pb, the muon binding energy is 
£6 = 10.475 MeV and the momentum transferred to 
the nucleus by the outgoing e~ of the /i^ — e~ con- 
version ranges from q = 0.482 fm^^, when the tran- 
sition energy E vanishes (namely in the coherent pro- 
cess), to zero, when E reaches the maximum value, i.e. 

— — Cfc = 95.183 MeV (namely when all the 
available energy of the bound /i~ goes to a nuclear exci- 
tation). The particle threshold energy i?thr is 8.09 MeV 
in the case of the SkM* force. We have obtained results 
for L = 0,1,..., 6 and for natural parity, (—1)'^. The 
most important contributions to the incoherent transi- 
tion rate are expected from L < 4 T^]. In all cases we 
have used F = 0.2 MeV and r„iax = 17 fm. 



A. Incoherent transition-strength distributions 
vesrus excitation energy 

In Figs, irm the distribution RaL{E) is plotted as a 
function of E, for L = 0,1,2,4, and for natural-parity 
transitions. In the monopole case, L = 0, the Isoscalar 
(IS) Giant Monopole Resonance (GMR) is the main peak. 
For 7 and Z, there is considerable contribution com- 
ing from higher energies (20-35 MeV), i.e., the isovec- 
tor (IV) GMR region. For L = 1, the IV Giant Dipole 
Resonance (GDR) corresponds to the strength clustered 
around E k, 12 MeV. For W and Z, important contribu- 
tion seems to come from higher energies (above 20 MeV) , 
in particular, the IS GDR. For W exchange, the region 
below 10 MeV contributes significantly. In this region 
we find the oscillation of the neutron skin against the 
nuclear core (pygmy dipole resonance) |T5l |. The tran- 
sition density is isoscalar in the interior of the nucleus, 
while on the surface the proton contribution vanishes. In 
the quadrupole case, L = 2, the IS Giant Quadrupole 
Resonance (GQR) is the second peak, close to 11 MeV. 
The collective low-lying peak is strong as well. There is 
some contribution from energies higher than 15 MeV, i.e. 
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FIG. 1: The distribution RaL{E) in ^°®Pb for L = 0. Skyrme 
parameterization SkM* has been used. 
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FIG. 2: Same as Fig. [T] for L = 1. 



from the IV GQR region, especially in the cases 7 and 
Z. For L = 3 (not shown) the strength is mostly con- 
centrated in the collective octupole state at low energy. 
For i > 3, as shown in Fig. ^for L = 4, the calculated 
strength is quite fragmented and most of it lies below 
20 MeV. 

A finite momentum transfer qf can result in the exci- 
tation of overtones of giant resonances [23 . l25l | . With the 
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FIG. 3: Same as Fig[T] for L = 2. 
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FIG. 4: Same as Fig[T] for L = 4. 



exception of the IS GDR, such states lie typically above 
30 MeV, for ^"^Pb, where the corresponding momentum 
transfer is less than 0.33 fm~^. As a result, the possible 
contribution of such states was too small to be identified 
and visible on graphs. 

The dipole results presented in Fig. |2 were obtained 
with corrected dipole operators - see Sec. lIIIBl for more 
details regarding this special case. A small amount of 
spurious strength remains close to zero energy (labeled 




FIG. 5: Fraction of the total strength SaL, coming from states 
below the particle threshold vs. the multipolarity L. Skyrme 
parameterization SkM* has been used. (Lines are drawn to 
guide the eye.) 



"sp." on the figure), but it is well separated from the rest 
of the distributions. 

In Fig.[5lwe plot the fraction of the total strength SaL 
coming from states below the particle threshold (S'aL.thr) 
and in Fig. the fraction coming from states below 
20 MeV (<S'aL,20Mcv), vs. the multipolarity L. In the 
dipole case, the remaining strength of the spurious state 
is not taken into account when evaluating these fractions. 

We see that for low multipoles L = 0, 1, 2 only a small 
portion of the strength originates from energies below 
particle threshold. The trend followed is similar for all 
mechanisms and, as we have verified, independent of the 
interaction used. For even multipoles L = 2, 4 a big por- 
tion of the contribution is pushed to higher energies as 
compared to the neighboring odd ones. Another interest- 
ing feature is the fact that, for some multipoles (L = 
for photonic mechanism, i = 1 for W^-boson exchange), 
a significant portion of the strength comes from above 
E = 20MeV. 

We have also calculated the fraction of the total 
strength SaL, coming from states below 50 MeV, for L 
up to 6. In all cases, the fraction is practically equal to 
unity. This means that the discretized versions of RPA 
and QRPA for the examined nucleus ^°*Pb are safe to use 
if the energy cutoff is large enough to sufficiently account 
for transitions below this value. 
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FIG. 6: As in Fig. |K| fraction of strength SaL, coming from 
states below 20 MeV. 



B. Dipole strength and the spurious CM motion 

It is well known [13. IT^ that the 1~ excitations contain 
admixtures of the spurious excitation of the center of 
mass (CM) of the nucleus 



R 



1 ^ 



(13) 



corresponding to a situation in which the unexcited nu- 
cleus moves as a whole around the localized fictitious 
potential well. Normally, these spurious components are 
separated out by the RPA methods. However, the use 
of a truncated model space and non-self-consistent single 
particle energies in ordinary RPA and the other versions 
of QRPA destroys the translational invariance and inserts 
spurious excitations into the spectrum. Thus, the spuri- 
ous CM state is not completely separated from the real 
(intrinsic) nuclear excitations, and in addition its energy 
eigenvalue is not zero. 

In Continuum-RPA models with Skyrme interactions 
it has been possible to achieve a high degree of self- 
consistency, i.e. the same interaction is used for the HF 
calculation of ground state properties and for the residual 
interaction. In addition, no truncation is involved. How- 
ever, due to the formulation of the model in coordinate 
space, it is common practice to exclude the Coulomb and 
spin-orbit contribution (at least) to the residual interac- 
tion. Therefore, self-consistency is violated and, even in 
cases where the spurious state appears very close to zero 
energy, some spurious strength may remain at higher en- 
ergies. 



For electric dipole excitations, the problem is usually 
treated by using effective charges ||2^. Similarly, in the 
case of IS dipole excitations, effective operators are used 
|l7l , which minimize the spurious admixture in the 
strength distribution. The effect on the IS dipole excita- 
tions of ^"*Pb was examined in detail in Ref. 28]. Here 
we will present a similar prescription for the operators 
involved in fi~ — e" conversion. 

We begin with a dipole excitation operator of the 
generic form 



= E^(^J-)/('^.-mo(r,) = ^c(r,) 




(14) 

As far as intrinsic excitations are concerned, this operator 
is equivalent to the corresponding "corrected" operator 



where stands for the z— component of the CM space 
vector R of Eq. (O and 




Our task is to determine the parameter rj so as to 
eliminate the spurious CM excitation. This can be 
achieved within the collective model, by imposing the 
translational-invariance condition on the transition den- 
sity characterising the collective state induced by ^f" 
]VA l2q . A condition on rj is thus obtained analytically. 
The result is (cp = c(l/2), c„ = c(-l/2)): 



1 



A dr 

where the mean values are defined as 



(16) 



(£'(?-))p,n = 



/o°°Pp,n(0^^dr 



and Pp{r), Pnir), are the proton, neutron, densities in 
the nuclear ground state. They are normalized so that 
4^ C PpA^y^^r ^ Z, N. 

By recalling from Eq. Q the operator Tai which in- 
duces the dipole 1~ excitations in fi~ — e~ , we find that 
it can be cast in the form 114|l with 

/(r) = 6gvfvVT2TTji{qr) 



and Cp = Ca(l/2) = Cap, c„ = Ca(-l/2) = Can- Thcu, 
using recursion relations of the Bessel functions, we find 
that 



J__d 

j,2 f^j. 



[f{r)r^] = 6gvfvVT2Trqjo{qr). 



(17) 
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Notice that the proton, neutron form factor Fp^n{q) = 
{jQ{qr))p,n- Under the above chxumstances Eq. 116|l 
reads 



A 



A 



-Fn{q) 



(18) 



The form factors are calculated numerically using the HF 
ground-state densities. 

The above prescription is a generalization of the 
method used in Ref. where a purely isoscalar field 
(cp = c„) was assumed and a specific form of the function 



/(r) was utihzed, i.e. /(r) 



In Eq. Hlt)|l the values of 



Cp,n and the form of /(r) are arbitrary. However, for the 
above mentioned isoscalar field, the present prescription 
leads to that given in Refs. [1II23. 

In Fig. [7| we plot the dipole distributions Rai{E) of 
Eq. ((SI, for photon, W- and Z-boson exchange diagrams, 
calculated by using the corrected and uncorrected oper- 
ator (they have been obtained with the SkM* force and 
for F = 0.2 MeV). One can see that most of the spurious 
strength below «6 MeV has been removed. The strength 
distributions above 20 MeV are practically unaffected. 
The strength between 6 and 20 MeV appeares some- 
what redistributed. The effect of the correction appears 
strongest in the case of the W^-boson exchange mecha- 
nism. We should note that the radial mesh used in the 
CRPA calculation, with Ar = 0.34 fm, may not be fine 
enough to yield completely converged results in this en- 
ergy region p^ . Numerical inaccuracies of such origin 
may be the reason why strength appears to be added at 
around 13 MeV (M^-boson case), rather than removed, 
after using corrected operators. A result that persists 
when other interactions are employed, is that the pygmy 
dipole state below 10 MeV is strongly affected by the 
correction in the photonic and W cases. 

In Table we list the portion of transition strength 
removed from the total contaminated 1^ transition 
strength Sai (stot)i as well as the portion of the strength 
removed from above 6 MeV excitation energy s^gM^v 
(with respect to the uncorrected strength above 6 MeV), 
for SkM* and for the three examined channels. For all 
three mechanisms, about 90% of the total transition rate 
was spurious. We expect this result to be independent of 
the interaction used, because the spurious state at low 
energy always dominates the isoscalar dipole strength 
distribution (which contributes in all three mechanisms) 
and because the corrected operators are, by construction, 
most effective for this state (thus removing practically all 
its strength). We were not able to demonstrate this in the 
particular cases of SGII and MSkT, because the energy 
of the spurious state in these cases was found imaginary, 
within the present calculation. In other words, we were 
not able to evaluate and take into account properly the 
strength of the spurious state, before or after the correc- 
tion. 



From Table n we notice that for s^bmcV ^'^^ 




the result is small in absolute value, but negative. It rep- 
resents the numerical accuracy of our calculation and, be- 
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FIG. 7: The dipole distributions Rai{E), for 7-photon and 
W^-boson exchange diagrams of the fi~ — » e~ conversion in 
^°*Pb. The results have been calculated for dipole opera- 
tor Tai (dotted line). In order to estimate the spurious CM 
contribution of this operator, we also show the distribution 
strength (full line) of the corresponding corrected operator 
given by Eq. HSU . 



ing small, it indicates that the degree of self consistency 
reached by our HF-I-CRPA model is sufficient to achieve a 
satisfactory separation of the spurious transition. For the 
W and Z cases, however, spurious admixtures of more 
than 6% are found above 6 MeV. These numbers vary 
when different Skyrme interactions are used, with their 
values remaining below 10%. (As mentioned before, they 
are not free of numerical inaccuracies.) One should ap- 
ply the same treatment in the case of other, not self- 
consistent RPA methods, where the lowest 1~ state is 
shifted artificially, by means of additional parameters, to 
zero energy. It is possible that larger corrections would 
be obtained, above 6 MeV. Such a result would mean 
that excluding the lowest 1~ state from the calculation 
of the incoherent rate, as was done on Ref. would 
not be an adequate treatment. As can be seen from Fig. 2 
of Ref. ^3 1 a considerable amount of spurious strength 
is distributed at higher excitation energies in the case of 
the QRPA calculations reported there. 



IV. SUMMARY AND CONCLUSIONS 

In the present paper we have focused on the investi- 
gation of the incoherent rate of the exotic /i^ — con- 
version in the heavy nuclear target ^°*Pb. We employed, 
for the first time in this process, the Continuum-RPA 



8 



fj. — e mechanism 


7 


W 


Z 




86.9 


96.3 


90.5 


®>6 Mcv(^) 


-1.3 


7.8 


6.1 



TABLE I: Percentage of the total 1^ transition strength Sai 
(s ^ot) and of the strength above 6 MeV mcv) consumed 
by spurious transitions, for the interaction SkM*, and for the 
three channels 7, W , Z. 

method which is appropriate for exphcit construction of 
the excited states lying in the continuum spectrum of the 
nuclear target. We mainly focused on the distribution of 
the transition strength as a function of the excitation 
energy of the target for energies up to 30 MeV. 

We have investigated in detail the various transition 
strength distributions coming from natural-parity ph ex- 
citations up to 1/ = 4 by using two different Skyrme in- 
teractions. We found that a significant portion of the 
incoherent fi~ — e~ rate comes from high-lying nuclear 



excitations. A similar study could be done for unnatural- 
parity transitions. 

The spurious 1^ admixture was eliminated by con- 
structing the purified dipole operators of the fj,~ — e~ 
conversion within the collective model. This study pro- 
vided us with the interesting result that the greatest por- 
tion of the 1~ transition strength is due to the spurious 
CM excitation, a result in agreement with that of an ex- 
act method constructed recently [T^ for removing spu- 
rious contaminations. The latter is a significant result 
for the fi~ — e~ conversion experiments searching for the 
coherent rate. 
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